Dynamical Anomalies and Intermittency in Burgers Turbulence 
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We analyze the field theory of fully developed Burgers turbulence. Its key elements are shock fields, 
which characterize the singularity statistics of the velocity field. The shock fields enter an operator 
product expansion describing intermittency. The latter is found to be constrained by dynamical 
anomalies expressing finite dissipation in the inviscid limit. The link between dynamical anomalies 
and intermittency is argued to be important in a wider context of turbulence. 
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A field theory of hydrodynamic turbulence is difficult 
in two ways: It is far from equilibrium and far from the 
realm of standard perturbative renormalization Q . New 
non-perturbative concepts have been established for sim- 
pler model systems sharing some characteristics of turbu- 
lent fluids. In particular, the stochastic Burgers equation 



d t \ + (v • V)v = v V 2 v + V?y 



(1) 



has become an important model for turbulence, recently 
studied by a variety of methods P~pl|. It governs the 
time evolution of a vortex-free velocity field v(r, t) = 
Vh(r,t). The equivalent scalar equation 



(2) 



d t h = is\7 2 h + -(yh) 2 + i] 



is known as the Kardar-Parisi-Zhang equation |12| , |13| . 
The driving potential r/(r,t) is Gauss distributed with 
mean rj(r, t) = and correlations 



r)(r, t)r)(r>,t>) = er 2 S(t - t')G(\r - r'\/r ) 



(3) 



over a characteristic scale ro . The function G is taken to 
be analytic with G(p) = 1 - G 2 p 2 + 0(p A ) and G 2 > 0. 
There is a second characteristic scale, the dissipation 
length ao = zy 3 / 4 ^ 1 / 4 . Burgers turbulence occurs at 
high values of the Reynolds number R = (r /a ) 4 ^ 3 and 
shows strong intermittency. For example, the longitudi- 
nal velocity difference moments of order k = 2, 3, . . . take 
the form 



(h(ri)-«||(r a )]' 



|n -r 2 \- kx - +£ To £k (4) 

in the inertial scaling regime tq/R <C |ri — r 2 1 <C r a 
The third moment grows linearly with |ri — r2 1 , i.e., 
x v = —1/3 and X3 = 0. The other moments, however, 
acquire a singular dependence on ro, which defines the 
intermittency exponents Xk- The turbulent state is asso- 
ciated with a particular strong-coupling limit of Eq. (Q) , 
called the turbulent limit in the sequel: v — > at fixed 
driving given by (|^). In that limit, the velocity field ac- 
quires discontinuities called shocks (see Fig. 1). It is these 
singularities that cause intermittency. 



Burgers' equation has a number of further applications. 
It has been proposed as a model for galaxy formation jl4| . 
The Kardar-Parisi-Zhang equation (H) models stochastic 
surface growth described by the "height" field h(r, t), and 
is related to directed polymers in a quenched random 
medium. In these contexts, the driving force is usually 
taken to be random in time and space, 

?7(r, t)r]{r',t') = j 2 5(t - t')5{r - r') . (5) 

For is 3 /"/ 2 — > 0, there is again a strong coupling limit, 
which, however, is quite different from the turbulent 
limit. The scaling is probably non-intermittent, 



((Mri)-Mr 2 ))' £ }^|r 1 -r 2 | 



— kx h 



(6) 



The exponent Xh depends on the dimension d. The value 
Xh = —1/2 for d = 1 is well-known |l5| ], while Xh = —2/5 
for d — 2 has been obtained only recently [ p^l3] | . 

This Letter studies the field theory of stationary Burg- 
ers turbulence. The basic quantities of this field theory 
are scaling fields containing powers of the local veloc- 
ity v and its gradients, Vv, VVv, etc. In particular, the 
statistics of the shocks is represented by a family of renor- 
malized shock fields (k = 1,2,.. .) that remain finite in 
the turbulent limit. Another family of fields Sk describes 
velocity gradients away from shocks. The short-distance 
properties of the scaling fields are encoded in an operator 
product expansion (OPE). This is a familiar concept in 
field theory (see, e.g., ref. |17|). It has recently been 




FIG. 1. Stationary Burgers turbulence in one dimension. 
The velocity profile v(r,t) at a given time t is an ensemble of 
shocks with amplitude ~ ito, width ~ a, and mutual distance 
~ ro. The probability distribution P r (u) of the longitudinal 
velocity difference over a distance a r -C ro is generated by 
no-shock and single-shock configurations (lines with arrows). 
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extended to non-equilibrium systems. OPEs have been 
proposed for general turbulent systems by Adzhemyan, 
Antonov, and Vasil'ev [18|, Eyink [|l9| (see also Du- 
plantier and Ludwig |^C|), and by L'vov and Procac- 



cia 1 21 , for Burgers turbulence by Polyakov |], and for 
the Kardar-Parisi-Zhang equation by Lassig |lq|. The 
Burgers OPE discussed here differs from Polyakov's con- 
jecture by the explicit inclusion of the fields and Sk- 
It is an important tool to understand the dynamics of 
Burgers turbulence. The shock singularities preventing a 
straightforward evaluation of the equation of motion, the 
stationary state is maintained by a rather subtle balance 
of driving forces, convection, and dissipation. We dis- 
cuss the resulting distribution of velocity differences (the 
functional form of which has been much debated p PI), 
as well as the flux of the energy density and its general- 
izations. The energy dissipated per unit of volume and 
time remains finite in the turbulent limit v — > 0. In 
the field theory of Burgers' equation, this is reflected by 
dynamical anomalies, i.e., asymptotically finite dissipa- 
tion terms in effective equations of motion for the inertial 
regime. These are associated with conservation laws that 
are valid for the inviscid equation without driving (y = 0, 
e = 0) but are broken in the driven state at any finite 
v |2^,^1. We find anomalies given by operator products 
involving the shock fields. This is not surprising since 
dissipation takes place at the shocks. Dynamical anoma- 
lies and intermittcncy are indeed closely related: they 
are both generated by the singularities of the velocity 
field. This conceptual link and the underlying theoreti- 
cal framework are expected to extend to other turbulent 
systems, as we briefly discuss at the end of this Letter. 

The basic phenomenology of stationary Burgers tur- 
bulence is well established in one dimension |3||l(J but 
appears to be much the same in higher dimensions . 
The velocity profile v(r, t) at a given time t looks similar 
to that of decaying Burgers turbulence with random ini- 
tial data |2^ |2jJ. It consists of ramps (i.e., regions where 
the velocity derivatives Vv, VVv, etc. are finite) sepa- 
rated by shock singularities, as sketched in Fig. 1. Shocks 
develop out of ramp regions through preshock singular- 
ities (which have the cubic root form v(r) — v(r ps ) ~ 
— (e(r — rps)) 1 / 3 for d = 1 H). The shocks have ampli- 
tudes of order Uq = {sTo) and distances of order tq 
to their neighbors; the slopes of the ramps are of order 
uq/to- These scales are independent of v, while the typ- 
ical shock width a = Tq/R = v/uq vanishes in the turbu- 
lent limit. From these shock characteristics, the velocity 
statistics can be inferred in an approximate way that ap- 
pears to give the correct scaling. Consider, for example, 
the longitudinal velocity differences = uii(ri,t)— uii(r2,i) 
(with vn = v • (ri — V2)/r and r = |i"i — r2 1) and 
the local excess velocity w(r,i) = v(r, t) — (v). These 
have normalized probability distributions P r (u)du and 
Q(w)d d w, respectively. Q(w) depends only on the ab- 
solute value |w| by rotational invariance and parity but 



P r (u) is an asymmetric function of u since the dynamics 
is not invariant under the transformation v — > —v. Both 
distributions are invariant under Galilei transformations 
v(r, t) — > v(r — vo<, t) — vo. In the turbulent limit, they 
are independent of v, i.e., of the shock width a. Hence, 
they can be written in scaling form, 

Pr(u) = ±7>(±,Pl, QW = \q(&), (7) 



u r r 



where u r ~ uor/ro is the average velocity increment over 
a distance r on a ramp. By Eq. (]?]), the single-point 
moments scale in a simple way, 



Xk 



-jfe/3 



(8) 



For p = r/r <C 1, the powers of the velocity difference 
can be expanded in the number of shocks present between 
ri and r 2 . One obtains u k ~ p Q u k + PiUq + 0(p2) with 
the n-shock probabilities po ~ 1 — p, p\ ~ p + 0{p 2 ), 
p n = 0(p n ) . This implies the bifractal moments 



(u k ) 



P t r -a/3 



(fc-3)/3 



(k<l) 
(fc>l) 



(9) 



What does all this mean for the field theory of the 
turbulent state? Correlation functions in the inertial 
regime should be represented by a field theory with short- 
distance cutoff a. Scale invariance emerges in the turbu- 
lent continuum limit a — > 0. Since that limit is nonsingu- 
lar for velocity correlations, the scaling dimension of the 
field v(r, t) takes the Kolmogorov value x v = —1/3; the 
fields \r k (r,t) have dimensions Xk = kx v (k = 1,2,.. .). 
Indeed, the distributions (Q) and the resulting moments 
(||), (0) are covariant under the scale transformations 
r — » br, tq — » bro, v — > b Xv v at fixed e. Turbulence does 
generate short-distance singularities for the moments of 
velocity gradients. Defining the longitudinal gradient 
V • v = u/r with a discretization length r ~ a, we write 



(V-v) fe (M) 



s*(r,t) 



(10) 



The fields Efe(r,t) and Sfe(r,t) represent the contribu- 
tions from configurations with and without a shock in 
the discretization interval. Using Eq. ([)]), we have 



,k-l 



Xk 



(sfc)rog ~ {uo/r Q ) k - r ' 



fc-3 
2fc/3 . 



(11) 
(12) 



The multiplicative renormalization factors Zk ~ a ab- 
sorb the short-distance singularities in ( ]Tl| ) , which defines 
the cutoff- independent fields of dimension ik- The 
slope fields Sk have different dimensions x' k determined 
by the regular part ( p^ ) of the ramp slope moments. 

We now make the assumption that the velocity field 
satisfies an OPE of the form 



2 



v(ri) . . . v(r fe ) 



-kx v -\-Xk 



C fc E fc (n) + 



(13) 



with r = |ri — r2|. Both sides of this relation are ten- 
sors of rank k whose indices are suppressed. The right 
hand side is a sum over all local scaling fields. The 
displayed terms contain the lowest-dimensional Galilei- 
invariant fields representing shock and ramp configura- 
tions, respectively, multiplied by dimcnsionlcss coefficient 
functions C k (r 12 /r, . . . ,r lk /r) and c k (r 12 /r, . . . ,r lk /r) 
(simple numbers for k = 2), and powers of r as required 
by dimensional analysis. The suppressed terms involve 
subleading Galilei-invariant fields and noninvariant fields 
such as v fc , Vv fc , etc. By differentiating ([l3|), one obtains 
a manifestly Galilei-invariant OPE for velocity gradients, 

k 

V • v(n) . . . V • v(r fc ) ~ \{ 8{v t - r i _ 1 )E fc (r 1 ) + 

i=2 

s k (rx) + ...; (14) 

see also [[l9| for general turbulent systems. Here the shock 
fields generate only contact singularities, while the slope 
fields have regular coefficients. It is not known whether 
there are other singular terms. 

The OPE is a consistency condition between 
renormalized correlation functions in the stationary 
state. It relates, for example, the n-point function 
(v(ri) . . . v(r n )) to the (n — k + l)-point functions 
(Efe(ri)v(r fe+1 )---v(r n )) and (s fe (ri)v(r fe+1 ) . . . v(r n )) 
in the limit a <C |r.y| <C |r^|,ro (i,j — l,...,k and 
I = fc+1, . . . , n). The coefficient functions C k etc. are as- 
sumed to describe local properties of the inertial regime 
independent of the cutoffs a and r$. Hence, the exis- 
tence of an OPE embodies two important characteristics 
of the turbulent field theory: (i) It is renormalized, i.e., 
the short-distance singularities of all scaling fields have 
been removed in a consistent way. (ii) The large-distance 
singularities are created solely by the single-point ampli- 



tudes of negative- dimensional fields, (O) ~ r 



In par- 



ticular, the intermittency exponent x k in Eq. (^J) equals 
the scaling dimension of the leading Galilei-invariant 
field in the expansion (13). This is indeed the case for 
Burgers turbulence, where (u k ) ~ (E^) ~ rQ Xk with 
ifc = — (k — 3)/3 according to Eqs. (||) and (pi]). 



The OPE ( J13[ ) has important consequences for Burgers 
dynamics, which we now discuss for simplicity in d = 1. 
For purely convective dynamics (y = 0, e = 0), the 
moments of the excess velocity w k (r,t) are locally con- 
served, d t w k + (k/(k+l))V(w k+1 ) = (fc = 1,2,. ..). In 
the driven state, (w k (r,t)) is pumped with a finite rate 



e(dlw k (r,t)) ~ r Q 



(k = 2,4,...). This cannot be 



offset by convection since (V(w k+1 )(r, t)) = 0. Hence, 
the stationary state must be maintained by a dynamical 
anomaly 



It is easy to check that the OPE ([L3|) predicts the correct 
form of the anomaly. We have 

^((V^Xr)^- 1 ^)) ~i/<*(0) (* fc (r)> 



using v 



4/3 



,V3 



e(d 2 w w k (r)) + v((V 2 w)(r) d w w k (v)) = 



(15) 



and a regularization on the scale 
of the short-distance cutoff, 6(0) ~ a -1 . Hence, the OPE 
and the resulting anomalies provide a link between the 
fields v k and the shock fields T, k . This fixes the intermit- 
tency exponents by the scaling relations x k — x k -2 + 1/3 
in accordance with (^) and ([ll]) . 

The role of dissipation is more subtle for the velocity 
differences. Consider again the shock number expansion 
P r {u) = (1 - p)P°(u) + P P^(u) + 0(p 2 ) for p < 1. By 
virtue of the OPE (|l3|), the equation of motion for the 
conditional distributions P^{u) and P^iu) is then essen- 
tially reduced to that of the families of single-point am- 
plitudes (E/c) and (s k ) Teg , respectively. 

1. The zero-shock part P^{u) is also covariant under 
the 'convective' scale transformations r — * br, u — > bu at 
fixed e and ro and can hence be written in the scaling 
form Pr(u) = (1 / u r )P° (u / u r ) discussed in ||. In other 
words, the expansion u k = r k s k + (l/2)r k+1 \7s k + ... is 
dominated by the first term, which depends on the large- 
distance scale ro only through the ramp slope moment 
(sfe) ~ (uo/ro) k . The convective symmetry is expected 
to be broken for —u^u r = (er) 1 / 3 by curvature effects 
at preshocks. This is precisely the scale where P r (u) 
becomes dominated by the shock part; see ( p7| ) below. 
In the equation of motion for V° , dissipation can be ne- 
glected at all points of finite \7 2 v (and even at preshocks). 
Furthermore, driving and convection are local processes 
in velocity space. It is straightforward to show that they 
are represented by the differential operators G2d 2 V°(ui) 
and (O^lu 2 + cj)"P°(a;), respectively [jsTJlQl. In particu- 
lar, the term ujP°(ui) describes a change in measure due 
to convective squeezing or stretching of the ramps [flo|| . 
We will be interested in solutions with a positive average 
ramp slope, i.e., with a net gain in measure, (lu) > 0. The 
formation of shocks, on the other hand, produces a mea- 
sure loss L(oj) < 0. In the stationary state, the net loss 
offsets the convective gain, (lu) + (L(oj)) = 0. For con- 
sistency with known properties of P r (u) , the equation of 
motion (G2d 2 +d LO uj 2 +lu+L)V° = must have a normal- 
ized, positive solution which behaves asymptotically as 
V° - wexp(-w 3 /3G 2 ) for w > 1 |,| and V° ~P\~ 3 ~ a 
for lu < -1, with a > 0, most likely a = 1/2 ||. This 
requires L(u) = o(u>) for to 1 and L(u>) ~ auj for 
u> <C — 1- The 'anomaly' L(lu) can be associated with 
ultraviolet-finite operator products Ei (r)s k (r) pq] . It 
is not clear, however, whether the functions L(u>) and, 
hence, P°(oj) are entirely determined by the OPE. 

2. The single-shock part P^ (u) is expected to have the 
scaling form P}(u) = (p / u )'P 1 {-u / u ) + 0{p 2 ). The 
leading term depends on r only through the shock prob- 



3 



V(u>,p) 



(17) 



ability pi ~ p in accordance with the OPE (|13|) and the 
moments ([)]). V x {a) is the scaled shock size distribu- 
tion function. The equation of motion for a single shock 
produces a driving term d 2 P 1 (a) and a convection term 
d a f (a)aP 1 (a) , while dissipation can again be neglected. 
Here (uo/ro)/(a) is the expectation value of the aver- 
age ramp slope to both sides of a shock of size uo<r. For 
very large shocks (<r ^> 1), this will be positive and pro- 
portional to the shock size, leading to the asymptotic 
equation of motion |^] (d 2 + AdaO^yP 1 = with a con- 
stant A > 0. This determines the tail of the shock size 
distribution, V 1 ~ exp(— A<r 3 /3) for a ^> 1, in agree- 
ment with instanton calculations j?J , while the dynamics 
of initial shock growth suggests V 1 ~ cr for a <C 1 0. 

The distribution P r (u) should be dominated by the 
zero-shock part for small velocity differences \u\ ~ u r 
and by the single-shock part for large negative values u ~ 
— uq. The crossover scale —u r = —u r p~ 2 / 3 is obtained 
by matching the two expressions, P®(—u r ) ~ pP^{—u r ). 
Thus we conjecture, using the scaling form of Eq. (Q), 

V°(u) (tV /3 > - 1) 

(?P l {-wp) (ujp 2 / 3 <-l) ■ 

It is easy to verify that this solution is indeed normaliz- 
able and has the correct moments (|9|). It has (u) > over 
ramps and (u) < at shocks, which is compatible with 
the constraint (u) = due to translational invariance. 

To summarize: Burgers field theory contains two differ- 
ent families of local scaling fields, and Sk, which rep- 
resent powers of singular and regular velocity gradients, 
respectively. The fields live on the shocks; their am- 
plitudes (S/t) generate intermittency. They are coupled 
to the other scaling fields through an OPE. The result- 
ing dynamical anomalies fix the intermittency exponents 
through scaling relations. The anomalies for velocity dif- 
ferences have a simple physical cause: the formation of 
singular velocity configurations out of regular ones. 

How much of this framework is preserved in Navier- 
Stokes turbulence can only be conjectured at present. 
Intermittency is still created by the infrared-divergent 
one-point amplitudes of negative-dimensional Galilei- 
invariant fields. These are associated with the singular- 
ities of the turbulent flow. Vortex filaments, for exam- 
ple, could play the role of the Burgers shocks [Q. The 
singularities have a much more complicated statistics, 
however. They lead to multifractal instead of bifractal 
scaling and may indeed suppress a coherent convective 
scaling regime P r (u) ~ (l/u r )P°(u>)- A stronger break- 
down of convective symmetry may produce anomalies 
and, hence, scaling relations compatible with intermit- 
tency exponents Xk nonlinear in k. Will such relations 
actually determine the values of Xk , leading to a nonper- 
turbative theory of turbulence? 

I thank Victor Yakhot for useful discussions. 
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